A Vertex Correction in the Gap Equation for the High Temperature Superconductors 
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We show that the Migdal theorem is obviously violated in the high Tc cuprates and the full vertex 
correction should be included, in particular, in the gap equation, in order to be consistent with 
the anomalously strong inelastic scattering in the region of the "hot spots" , which is observed from 
the various normal state experiments. The full vertex correction is obtained non-pertubatively by 
utilizing the generalized Ward identity, which is shown to hold in the important scattering channel 
of the pairing interaction in the high Tc cuprates. As a result, we find a strong enhancement of Tc 
from the vertex correction despite of the strong pair breaking effect due to the inelastic scattering. 
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Since the discovery of the high Tc superconductors, 
there has been substantial progress in our understand- 
ing for the large part of normal state anomalies, if 
not its microscopic origin. The marginal Fermi liq- 
uid (MFL) ansatz [Q, for example, captures remarkably 
simple essence to understand the various normal state 
anomalies in a unified manner. The single and essential 
ingredient of the MFL phenomenology is the assumption 
of the anomalous scattering kernel only with tempera- 
ture scale. Although the microscopic mechanism of this 
scattering is yet controversial, it is generally agreed on 
its existence and the essential role of it for the high tem- 
perature superconductivity (HTSC) phenomena. More 
recently, however, it becomes clear that the MFL ansatz 
is rather too simple to explain more problematic normal 
state anomalies, in particular, such as Hall resistance, 
c-axis charge dynamics, etc. The angle resolved photo 
emission (ARPES) experiment for example, reveals 
strong anisotropy of the scattering rate along the Fermi 
surface, conveniently termed as the "hot spot" region and 
the "cold spot" region according to strong or weak scat- 
tering rates in the different sections of the Fermi surface. 

While it is natural to think that such an anomalous 
inelastic scattering, which is responsible for the normal 
state anomalies, should also play a role in the supercon- 
ducting pairing mechanism, there has been no satisfac- 
tory attempt to incorporate it in the pairing mechanism 
beyond the leading order. Up to now, the superconduct- 
ing pairing mechanisms for the high Tc cuprates can be 
largely grouped into two: non-Fermi liquid theories and 
Fermi liquid based theories. The former is typically rep- 
resented by the interlayer pair tunneling mechanism || 
and the latter one has more variety such as the anti- 
ferromagnetic paramagnon theory Q, Van Hove singu- 
larity, or more exotic theory like various preformed boson 
scenarios ||. Each theory has various degree of success 
to account the known data and we are not in the posi- 
tion to judge of them. In this paper, we would like to 
address the importance of the vertex correction when we 
apply the conventional BCS-Eliashberg pairing approach 
for HTSC. 

It has been pointed out by many authors the inade- 
quacy of the blind application of the Migdal theorem |i]|7]] 
in HTSC. At the qualitative level, there arc already quite 



simple reasons not to trust the Migdal theorem in HTSC: 
(1) wo/Ep is not so small (wn is the characteristic fre- 
quency of any mediating boson); (2) the Fermi surface 
anisotropy (certain degree of nesting or Van Hove singu- 
larity) poses a potential danger to invalidate the phase 
space constraint in the Migdal theorem argument. Nev- 
ertheless, the complexity of higher order vertex diagrams 
prohibits systematic progress along this line; at most it 
can only indicate the danger of vertex correction by cal- 
culating the first order correction. 

In this paper, we took rather simple short cut by ob- 
serving two experimental findings: (1) HTSC is a d-wave 
pairing state mediated by strongly peaked potential in 
the momentum space, which dominantly mediates scat- 
tering from one "hot spot" region to the other "hot spot" 
region in the Brillion zone ||; (2) the electrons in the 
"hot spot" region has a singular self-energy correction, 
destroying almost its quasi-particle nature at all. From 
the above two observations, we show that we are exactly 
opposite limit, for the important pairing channel, from 
which the Migdal theorem is valid. In this opposite limit 
- we conveniently call it the "Ward identity limit" - we 
can easily read off the exact vertex from the generalized 
Ward identity || , given a set of reasonable assumptions. 
As a demonstration, we performed a numerical calcula- 
tion to solve the model Tc equation for a d-wave state 
with a full vertex correction. We find a strong enhance- 
ment of Tc from the singular vertex correction despite of 
the strong pair breaking effect of the self-energy correc- 
tion. 

We briefly reexamine the Migdal theorem. Fig.l.a 
shows the typical vertex correction diagram of the first 
order. The typically interesting parameter regime of the 
momenta is that q,p,k ~ kp and qo,po,ko ~ wd] here 
and afterward wd is not necessarily the Debye frequency, 
but some characteristic frequency of any mediating bo- 
son. For the Migdal limit, i.e. q '■ vp > qo(^ wrj), the 
correction is of the order 0{wd/Ef) either due to the 
energy denominator in the electron propagators or due 
to the phase space constraint. Here the important obser- 
vation is that q-VF should be understood as e(k+q) — e(k) 
in more general expression and particularly for the tight 
binding band. For the other limit, i.e. q ■ vf < 90j 
which we call the "Ward identity limit" , the correction 
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is O(l) as found by Migdal and also more clearly demon- 
strated by Engelsberg and Schrieffer || using the gen- 
eralized Ward identity. However, the above conclusion 
is true only under the following set of assumptions: (1) 
electron-phonon vertex, g(k,k'), has no strong momen- 
tum dependence - in practice, assumed to be constant 
in most of analysis; (2) similarly, the phonon dispersion, 
w q , is assumed to be isotropic; namely, it doesn't have 
any special directionality; (3) the Fermi surface doesn't 
have a strong anisotropy so that it doesn't introduce any 
special scattering phase space constraint; a pathological 
case is one dimensional system, for example; (4) in par- 
ticular, for superconducting instability, the typical scat- 
tering process is involved with large momentum exchange 
(\q\ ~ 0(|fc^|)); for a contrasting example, the dynamic 
polarizibility, P(q 7 w), of electron gas for \q\/pF — * limit 
should have an important vertex correction as a trivial 
violation of the Migdal theorem Q . 

Now let us examine the situation of the high Tc 
cuprates. First of all, wp/Ep is not that small since 
now the relevant wp is most probably electronic origin 
and also Ep is renormalizcd to a smaller value |9j. How- 
ever, it is just one general fact to invalidate the Migdal 
theorem and it is not our main concern in this paper. 
There are two important observations just specific to the 
high Tc cuprates which break the Migdal theorem. First, 
the hole-doped high Tc cuprates is known to be a d- 
wave superconductor by now and accordingly the pair- 
ing potential V(k, fc') has a strong momentum depen- 
dence; we have in our mind the anti-ferromagnetic spin 
fluctuation mediated potential which is strongly peaked 
for fc — k' = Q, where Q ~ (±7r, ±7r) Q. Second, the 
Fermi surface is anisotropic; in particular, the impor- 
tant parts of the Fermi surface with the maximum gap 
opened are connected by Q (see Fig.l.b). Combining 
these two facts violates the Migdal assumptions (l)-(3) 
above. In particular, the scattering phase space con- 
straint - the relative probability for satisfying q-vp < wd 
with \q\ ~ \kp\ is 0(wp / Ep) with the Migdal assump- 
tions - doesn't hold anymore. Moreover, the last condi- 
tion for superconductivity of large momentum exchange 
doesn't help, too, because the superficial Migdal condi- 
tion, i.e., q ■ vp > wp> with large q for a parabolic band 
becomes e(k + q) — e(k) < wp, for the high Tc supercon- 
ductors when both k + q and k belong to the hot spot 
region. To summarize, in the high Tc cuprates, because 
of the strongly momentum dependent pairing potential, 
V^(fc,fc'), inducing a d-wave pairing and the anisotropic 
Fermi surface from the tight binding nature, the impor- 
tant scattering process for the superconducting pairing is 
not in the Migdal regime (e(k + q) — e(k) > wp) but in 
the Ward identity regime (e(fc + q) — e(fc) < wp). In the 
following, we show that we can indeed extract the full 
vertex correction of this scattering channel (fc, fc + q G 
"hot spot" region) from the generalized Ward identity 




Engelsberg and Schrieffer [|| derived the following gen- 
eral identity from the particle number conservation con- 
dition. 

q T Q (k; k + q)- q- f (fc; k + q) = G'^k + q) - G~\k) 

(1) 

where fc, q stands for the four momenta, i.e., k = (ko, fc), 
and r° and T are scalar and vector vertices, respectively. 
G(k) = (fco — e (fc) — ^(fc)) -1 is the full electron green 
function with a self-energy. The above equation is ex- 
act for any k and q and the famous "Ward identity" 
is derived by taking q = and qg — > limit, so that 
r°(fc;fc) = 1 - dE(k)/dk . The other limit, q = and 
q — > 0, can be used to obtain the vector vertex r(fc; fc). 
Since we assume a d-wave pairing in the high Tc super- 
conductors and the maximum gap is formed around the 
"hot spots" , the important part of the pair potential is 
the process which scatters from one "hot spot" region 
to the other "hot spot" region. Consequently, our inter- 
est is the scalar vertex, r°(fc, ko;q ~ Q,qo < wp), k G 
"hot spot" region. In general, there is no simple solution 
for this unless we have a simple form of £(fc). Theo- 
retical (ll) and experimental || studies show that the 
self-energy correction around the "hot spots" is singular 
in frequency but its momentum dependence is negligible. 
Therefore, we can separate the scalar and vector vertices 
to a good degree of approximation as follows, for the in- 
teresting regime of q ~ Q, qo < wd and k £ "hot spot" 
region J12| . 

g r (fc,fc ;Q,eo)= (2) 

Qo - [(S(fc +Q,k + q ) - S(fc, fc ))], 

Q ■ f (fc, fc ; fc + Q, fc + qo) = e(k + Q) - e(fc), (3) 

where e(fc) is the bare band dispersion. As argued above, 
we assumed that E(fc + Q, fco) ~ S(fc, fco) when both fc + 
Q and fc belong to the hot spot region for the above 
separation. 

Now to be specific, we approximate the self-energy for 
the hot spot region as the marginal Fermi liquid type as 
follows. 

ImYi{w) = a\w\ for \w\ < w c , (4) 

where w c is the high energy cutoff, and a is about 0.6 
from experiment |13| but it can be treated as a param- 
eter for our purpose. In reality, the self-energy correc- 
tion can be even more singular in the hot spot region 
p"l| , or maybe just Fermi liquid type (S ~ w 2 ) for 
very low frequency limit. In any case, our main con- 
clusion doesn't change, i.e. large vertex correction from 
the anomalous scattering in the hot spot region. Here the 
marginal Fermi liquid type self-energy assumption is just 
for demonstration purpose; nevertheless we think it is 
still reasonable assumption in view of experiments. Now 



2 



the full self-energy in the Matsubara frequency is writ- 



ten as E(w„) 
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- . And the corresponding 



scalar vertex is written as 
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T°(k,w n ;q ~ Q, Q„) 
[S(w„ + il n ) - £(w„)]/ifi„, 



(5) 



where w n = ~KkBT{2n + 1) and Q n = nkBT(2n) are 
the fcrmionic and bosonic Matsubara frequencies, respec- 
tively. Since we are interested in 0„ < wd regime, we 
take the Q n — > limit from the above equation. Now the 
full vertex for q ~ Q, k £ "hot spot" region in the static 
limit is the following. 



r°K) = [i + -in^ 



ti'7, 



(6) 



This vertex shows the InT divergence when w n — > limit 
due to the singular self-energy 14 1. The first term of 
the Eq.(6) is nothing but Z n in Matsubara formalism; 
Z n is defined as iw n — £(u>„) = iw n Z n . We notice that 
the singular suppression of the effective coupling constant 
due to Z n as g e ff = g/Z n can be almost recovered now 
as g e ff — gT^/Z n by the corresponding vertex correction 
reminiscent of the original Ward identity in QED. With 
this vertex we solve the model Tc equation for a d-wave 
state. The Tc equation is written as 



Z nK 



(7) 



Some remarks are in order about the above Tc equa- 
tion. First of all, the above Tc equation is the static 
limit of the Eliashberg equation, so that it would be the 
BCS equation were it not for the self-energy correction, 
Z n , and the vertex correction, T°(w n ) of static limit. We 
could have solved a full Eliashberg equation with a given 
dynamic pair potential and it would be no problem to 
include the dynamic vertex correction, T°(w n ;£l n ), from 
Eq.(5). We think that qualitative results would not be 
changed. Second, the vertex correction we put in the Tc 
equation is valid only when both incoming momentum, 
k, and outgoing momentum, k' = k + Q + Ak, belong to 
the hot spot region. In order to simulate this hot spot 
scattering channel constraint, we introduced a function 
$(fc, k'), which is just 1 only when both k and k' belong to 
the hot spot region connected by Q; otherwise it replaces 
T°(w n ) with the bare vertex, r° = 1. It means that we 
do not consider any vertex correction for other scatter- 
ing channels other than this special scattering channel. 
Some justifications for this are: (1) for other scattering 
channels, simply the vertex correction should not be as 
singular as for the hot spot scattering channel; (2) even 



if there should be some vertex correction, its effect in the 
Tc equation is suppressed by the pair potential, V(k, k'), 
since it is peaked only for the hot spot scattering channel. 

For simplicity of numerical calculation, we assume 
V{k,k') = V\sin(<f>-(f>')\ and A(fc) = A max cos(2(t>) in 2- 
dimensional momentum space, a circular Fermi surface, 
and <f> is the angle along the Fermi surface. For all cal- 
culations, we take w c = Q.5eV, wd — O.ieV (wd enters 
as the BCS cutoff in the momentum summation) and the 
coupling constant parameter A = VN(0) is taken to be 
1.5 (N(0) is the density of states at the Fermi level). This 
choice of parameters is only for exemplary purpose. Also 
for the hot spot scattering channel constraint function, 
$(fc,A/), we define the the hot spot region by the angle 
8 hot a s indicated in Fig. lb; now $(0, <^')=1 only when <f), 
4> G "hot spot regions" , otherwise it will set T°(w n ) = 1. 

In Fig2, we show Tc as a function of the strength of the 
self-energy correction, a, for different values of Ohot- For 
example, Ohot — tt/2 is the case that the whole Fermi sur- 
face is treated as the hot spot region. As expected, the 
8 hot — 7r /2 case (open square) shows an extreme slow- 
down of the suppression rate of Tc because of the over- 
imposed vertex correction; the solid square indicates the 
result with no vertex correction at all but only with the 
self-energy correction, showing rapid suppression of Tc 
with increasing a. In between these two curves, the re- 
sults with 9 ] w t — 7r/4 and 7r/8 are shown and these should 
be the more realistic cases. In short, Fig. 2 shows the dra- 
matic effect of the vertex correction for Tc even with a 
very narrow area of the hot spot region (see 9hot — tt/8 
case). When considering a ~ 0.6 for YBCO (Tc ~ 90 K) 
[ PH , without the vertex correction at all (solid square), 
Tc ~ 24 K only; the set of parameters (the dimensionless 
coupling constant, A = 1.5, and a characteristic energy 
of mediating boson, wd = 0-SeV ) is already quite favor- 
able choice for the pairing. Only with Ohot — 7r/8, Tc is 
enhanced more than by a factor of 3 (Tc ~ 85 K) for the 
same a = 0.6. The message is that in order to achieve Tc 
~ 100 K of a d-wave state in HTSC, with strong inelastic 
scattering but without including the corresponding ver- 
tex correction, we are forced to choose rather unrealistic 
parameters in the conventional pairing model. 

In conclusion, we show that the Migdal theorem is max- 
imally violated in HTSC, not only because of wd/ Ep ~ 
0(1) but more importantly because of the strongly mo- 
mentum dependent pairing potential, V(k,k'), and the 
existence of hot spot region in the anisotropic Fermi sur- 
face, which are dominantly scattered by the pairing po- 
tential in the main pairing channel of a d-wave state. 
Then, we show that we are in the "Ward identity limit" 
which allows us to extract the full vertex correction 
from the self-energy. Assuming a phenomenological self- 
energy (X = a|u>|) and taking into account the cor- 
responding vertex correction around the "hot spots" , we 
solved the model Tc equation for a d-wave pairing includ- 
ing both the self-energy correction and the vertex correc- 
tion. The results show the dramatic enhancing effect 
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of Tc by the vertex correction despite of the strong pair 
breaking self-energy correction. Considering the problem 
of choosing realistic values of parameters for a d-wave 
pairing model, our result enforces the essential role of 
the singular vertex correction, at least for a conventional 
BCS-Eliashberg type pairing scenario. Amusing obser- 
vation is that the strong inelastic scattering in the high 
Tc cuprates, which is the key entity responsible for the 
various anomalous normal state properties Q , turns out 
to be not so much destructive for superconductivity in 
the end, in contrast to the conventional superconductiv- 
ity, thanks to the singular vertex correction. The origin 
of it is the strong anisotropics in the pairing potential 
V(k, k'), the Fermi surface, and the order parameters 
A(k), which are all related, though. Therefore, in HTSC 
the strong correlations not only in frequency domain but 
also in spatial domain play important roles on their own 
to make all these unusual materials. 
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The hot spot region is defined as Afc£ < 1. Afc is the 
distance of k from the nearest "hot spot" along the Fermi 
surface. In optimally doped YBCO sample, £ ~ 2 above 
Tc, therefore the hot spot region can have substantial 
fraction along the FS (see, e.g., B. P. Stojkovic and D. 
Pine in Ref.[ll]) 

Z. Schlesinger et al., Phys. Rev. Lett. 65, 801 (1990). 
If a Fermi liquid like self-energy, say, £ = aw with a ~ 
0(1/ Ef), the vertex correction is still O(l) as w n — > 0. 
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FIG. 1. (a) Vertex diagram of the first order; solid lines 
are electron propagators and dotted line is any bosonic prop- 
agator, k and q stand for four momenta, i.e., k = (ko,k). 
(b) Schematic Fermi surface of HTSC; the areas around zone 
corners, (tt,0), and its symmetry related points, are the 
"hot spot" region, which are connected by the momentum 

Q ~ (±7r,±7r). The angle Q^ot defines the hot spot region 
quantitatively for the numerical calculation. 
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FIG. 2. Tc vs a, the pair breaking parameter defined in 
E = a\w\; solid square is the one without vertex correction, 
and open upper triangle, open circle, and open square are the 
ones with more vertex correction (wider area of the hot spot 
region) in increasing order. 
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